On the applicability of the equations-of-motion technique for quantum dots 
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The equations-of-motion (EOM) hierarchy satisfied by the Green functions of a quantum dot 
embedded in an external mesoscopic network is considered within a high-order decoupling approxi- 
mation scheme. Exact analytic solutions of the resulting coupled integral equations are presented in 
several limits. In particular, it is found that at the particle-hole symmetric point the EOM Green 
function is temperature-independent due to a discontinuous change in the imaginary part of the 
interacting self-energy. However, this imaginary part obeys the Fermi liquid unitarity requirement 
away from this special point, at zero temperature. Results for the occupation numbers, the density 
of states and the local spin susceptibility are compared with exact Fermi liquid relations and the 
Bethe ansatz solution. The approximation is found to be very accurate far from the Kondo regime. 
In contrast, the description of the Kondo effect is valid on a qualitative level only. In particular, 
we find that the Friedel sum rule is considerably violated, up to 30%, and the spin susceptibility is 
underestimated. We show that the widely-used simplified version of the EOM method, which does 
not account fully for the correlations on the network, fails to produce the Kondo correlations even 
qualitatively. 
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I. INTRODUCTION 

Quantum dots embedded in mesoscopic structures are 
currently of great experimental and theoretical interest, 
because such systems allow for detailed and controlled 
studies of the effects of electronic correlations^* The the- 
oretical description of these systems is usually based on 
the Anderson model^ in which the electronic correla- 
tions are confined to few impurities that represent the 
quantum dots. Although a rich variety of techniques 
has been developed over the years to treat the Anderson 
model^ their applications to the quantum dot systems 
are not straightforward. In such systems one would like 
to be able to study dynamical properties (e.g., transport) 
as function of the impurity characteristics (which can be 
tuned experimentally) over a wide range of parameters. 
This is not easily accomplished by the Bethe ansatz so- 
lution, for example. 

The theoretical difficulty can be pinned down to the 
ability to derive a reliable, easy-to-handle, expression 
for the single-electron Green function on the quantum 
dot. Because the electronic interactions in the Anderson 
model take place solely on the dot, this Green function 
can be shown to determine, under certain conditions, the 
charge or the spin transmission through the quantum dot, 
the charge accumulated on it, etci 3 ^ Consequently, much 
effort has been devoted to finding faithful analytic ap- 
proximations for this object. Alternative treatments rely 
on numerical techniques, such as quantum Monte Carlo, 6 
or the numerical renormalization group (NRG) method. 7 
NRG in particular is considered to be capable of provid- 
ing accurate estimates of the Green function over a wide 
parameter range, although at the cost of running itera- 



tive diagonalizations for each parameter set, and a lim- 
ited resolution at high energies and high magnetic fields. 

A ubiquitous method to derive an analytical expres- 
sion for the Green functions is to use the equations-of- 
motion (EOM) A In the case of the single- impurity An- 
derson model, the EOM of the (single-particle) impu- 
rity Green function gives rise to an infinite hierarchy of 
EOM of higher-order Green functions. A well-known ap- 
proximation procedure is then to truncate this hierar- 
chy, thus producing certain thermal averages represent- 
ing various correlation functions. The latter need to 
be found self-consistently from the resulting closed set 
of equations. The level at which the EOM are trun- 
cated is chosen such that most of the interaction effects 
are capturedAiSiiiii^ This scheme has been applied to 
the original Anderson model a long time a g 0f 2ii2iiAiiiii4 
yielding approximate expressions for the resistivityj^iifli 
the spin susceptibility^2*A£ and the magneto transport * 4 
of dilute magnetic alloys. For temperatures above the 
Kondo temperature, Tk, these results agree with pertur- 
bation theory calculations, 4 Although several limitations 
of these self-consistent approaches are known (such as un- 
derestimation of the Kondo temperature or the absence 
of (T/Tk) 2 terms in the low-temperature expansion of 
the results), they still can form a basis for a qualitative 
analytic treatment of the Kondo effect. 

The first application of the EOM technique to a quan- 
tum dot system has been undertaken by Meir, Wingreen 
and Lee^ (MWL). Neglecting certain correlation func- 
tions, they have obtained a closed analytic expression 
for the dot Green function. Several subsequent studies 
have employed their scheme to describe various effects of 
the Kondo correlations in quantum dots in different set- 
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tings, e.g., ac response of a biased quantum dot^ non- 
equilibrium Andreev tunneling^ and coupling to mag- 
netic leadstf As we show below, the approximation of 
MWL fails in various aspects. Recent attempts to im- 
prove on that solution turned out to be not completely 
satisfactory: either requiring further approximations on 
top of the self-consistent truncation^ or leading to it- 
erative numerical solutions>i2i2Si2i It is therefore of in- 
terest to examine whether an exact solution of the self- 
consistency equations of the EOM method will improve 
on the previous results, and will be capable of producing 
a reliable approximation for the Green function that can 
be used in the analysis of quantum dot systems. 



In this paper we consider a quantum dot embedded 
in a general complex network, and obtain and solve the 
truncated EOM for its Green function. This EOM con- 
tains all the next order correlation functions, which we 
decouple and calculate self-consistently. We first derive 
(in Sec. [nj an integral equation for the Green function, 
allowing for arbitrary values of the on-site Coulomb inter- 
action, U. We then analyze (in Sec. IHIjl the properties of 
this EOM solution. We find that in the case of a particle- 
hole symmetric Hamiltonian the resulting Green function 
is temperature-independent (a point which is overlooked 
in previous treatments of quantum dots, see for exam- 
ple Ref. S). 

However, we show that this is a singular 
point in the parameter space. We also investigate the 
zero-temperature limit of the EOM solution, and show 
that it fails to satisfy the Friedel sum-rule. We then turn 
to the infinite U case (in Sec. IIV|I . and derive an ex- 
act analytical solution to the integral equation for the 
Green function. In Sec. El we use this solution to obtain 
the total occupation number on the dot, and compare 
it with the exact solution supplied by the Bethe ansatz. 
This comparison shows that the EOM solution is faithful 
outside the Kondo regime, but fails in the regime where 
Kondo correlations play a dominant role. In particu- 
lar, the Friedel sum-rule is violated by ~ 30%, invalidat- 
ing the assumptions made in previous studies^A*i2i2I We 
then derive the local spin-susceptibility on the dot and 
demonstrate that the full self-consistent solution of the 
EOM, as used in this work, is required in order to obtain 
quantitatively correct results. We also examine the local 
density of states at the Fermi level, and find that it shows 
the expected universal behavior as function of T/Tk, 
though the EOM Kondo temperature lacks a factor of 
two in the exponential dependence on the single-energy 
level on the dot. Finally, we examine the EOM tech- 
nique from another point of view (Sec. IVIfl : We expand 
the Green function derived from the EOM to second or- 
der in the dot-network coupling, and compare the results 
with those obtained from a straightforward perturbation 
theory iSSiS This comparison shows again the necessity 
to include in the EOM solution the full self-consistent 
calculation of all the correlations. A short summary in 
Sec. IVIll concludes the paper. 



II. THE GREEN FUNCTION ON THE DOT 

As mentioned above, several properties of the Ander- 
son model can be expressed in terms of the Green func- 
tion on the dot. Here we examine the determination of 
this function using the EOM method. Our discussion is 
limited to a single interacting impurity embedded in a 
general non-interacting network, for which the Hamilto- 
nian can be written in the form 

H = Hdot + "H ne t + Wnct-dot ■ (1) 

Here, the dot Hamiltonian is given by 
Hdot = y^(eo + (rh)n da + Un d ^n dl , n d<J = d\d a , (2) 

where d\ (d a ) is the creation (annihilation) operator of 
an electron of spin index a = ±1/2 on the dot, eo is 
the single-particle energy there, h is the Zeeman split- 
ting, and U denotes the Coulomb repulsion energy. The 
non-interacting network is described by the tight-binding 
Hamiltonian 

Hnct = y t En(T a na a nij ~ ^ t ^mn a rncr a niT ' (*0 
no- rana 

where (a ncr ) is the creation (annihilation) operator 
of an electron of spin index a on the nth site on the 
network, whose on-site energy is e na , and J mn = J* m 
are the hopping amplitudes on the network. Finally, the 
coupling between the dot and the network is given by 

Hnct-dot = - J n<r d l a na + H - c - ( 4 ) 

We have allowed for spin-dependent on-site energies on 
the network, as well as spin-dependent hopping ampli- 
tudes between the dot and the network. In this way, our 
model includes also the case of spin-polarized leads con- 
nected to a quantum dot (see for example, Refs. and 
l2ll) . The entire system is assumed to be at equilibrium 
with a reservoir held at temperature T (in energy units) 
and chemical potential /i = 0. 

Adopting the notations of Ref. HI, we write a general 
Green function in the form 

/+00 
Q(±t)([A(t);B] + )e^ ± ^ t dt, (5) 
-00 

where A and B are operators, O is the Heaviside function, 
and 77 — > + . The Green function on the dot is then 

G a (z) = 4>*> z = uj±ir). (6) 

In conjunction with the definition l|5Jl. a thermal average, 
(BA), is related to the corresponding Green function by 

(BA)=i£^f(z)^A;B^ z , /(*) = —L^ , (7) 
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where the contour C runs clockwise around the real axis. 

The EOM for the dot Green function is given in Ap- 
pendix El [see Eq. (|A2|) ] . As is shown there, that equa- 
tion includes a higher-order Green function, whose EOM 
gives rise to additional Green functions. The resulting 
infinite hierarchy of EOM is then truncated according 
to a scheme proposed originally by Mattis^ and subse- 
quently used in Refs. Hl(miBl3ll4ll8ll9l2fl2ll25l Each 
Green function of the type <A t BC; dt» in the EOM 
hierarchy is replaced by 

^BC; rft> => (At-BKC; rft> - (A^C)^B; rft> (8) 
I 



if at least two of the operators A, B and C are network 
operators a n<7 . Explicitly, the Green functions which are 
decoupled are <a ntT a%a m ^, 4>, <al l& d g a m(T ; 4», and 
<^a} n ^a mS d a \ d\~^> (with a — — a)*L Upon calculating 
the averages {a\ na a m , a ) and (d\a m(T ) using Eq. 0, the 
set of EOM is closed, and can be therefore solved. The 
details of this calculation are presented in Appendix lAl 
In particular, the resulting equation determining the dot 
Green function is 



q ( s = U(z) - (n d9 ) - Pgjzx) ~ Pg(z 2 ) . , 

a[ ' u(z) [z- e -ah- E ff (js)] + [P 9 (zx) + P*(z 2 )} E ff (z) - Q 9 (zi) + Q s {z 2 ) ' 1 ' 

where 

u{z) = U~ l [U - z + e + ah + E ff (z) + E ff (z x ) - E ff (afa)] , (10) 



and z\ = z — 2ah, z 2 = — z+2eq + U . The functions P and 
Q are given in terms of the non-interacting self-energy 
on the dot, E CT (z), brought about by its coupling to the 
network [namely, the self-energy of the non-interacting 
dot, see Eq. (|AT|) ] . and the dot Green function itself [see 
Eqs. (H>J], 

P<,(z)=$„[G[ , Q„(z) =5„[1 + SG] , (11) 
where the notation ^ az [g] stands for 

ZvM = f / f{w)g a { W ) ^ {w) ~^ [z) dw. (12) 
2ir J c z — w 

Equa tion ffl gen eralizes the result of Ref. (see also 
Rcfs. 18 19 2(121(1 for the case in which the interaction 
on the dot is finite, and the entire system is subject to an 
external magnetic field. Our generalization also corrects 
a few details in Lacroix's earlier treatment of finite U. 12 



III. PROPERTIES OF THE EOM 
APPROXIMATION AT FINITE U 

As is evident from Eq. J3J, the solution of the dot 
Green function within the EOM scheme cannot be easily 
obtained over the entire parameter range. However, there 
are certain limiting cases in which this Green function 
can be analyzed analytically. We examine those in the 
subsequent subsections. 

A. Particle-hole symmetry 

Upon replacing the particle operators by the hole ones, 
d a = d rT ,al ia = a na , the Anderson Hamiltonian Q at- 



tains its original structure, with 

?o + ah = —€q — ah — U , U = U , 

J na J'na ) ^na — ^ncr 3 Jjim — Jmn ■ (13) 

(Hole quantities are denoted by a tilde.) The dot Green 
function in terms of the hole operators is then related to 
the particle Green function by 

G a {z) = <J CT ; - -G a {-z) . (14) 

One may check that this equivalence holds by introducing 
the definitions ifT^jl into Eq. Since E CT (z) = —T, a (—z) 
and u(z) = 1 — u(—z), one finds that [see Eqs. (|A25|) 
and SKM ] P a (z) = -Pa(-z) and Q a (z) = Q a {-z) - 
S CT (—z), re-confirming Eq. i|14|) . 

From now on we shall assume that Yi a (z) = —Y, a (—z). 
This relation is realized, for example, when the network 
to which the dot is coupled has a wide band spectrum, 
with the Fermi level at the middle. We next discuss the 
particular point where 2eo + U = and h = 0. At this 
point, the Anderson Hamiltonian becomes particle-hole 
symmetric. Then G„(z) = —G a {—z) 1 P<r(z) + P a (—z) = 
0, Qa{z) - Qa(-z) = £ CT (z), and (n d<J ) = 1/2. As a 
result, Eq. © becomes 

[GA-AY 1 - 



4[z-E ff (z)-2S ff (z)] ' v ; 

Namely, at the particle-hole symmetry point of the 
Anderson model, the EOM results in a temperature- 
independent dot Green function! This implies that the 
EOM technique at the particle-hole symmetric point can- 
not produce the Kondo singularity. This property of the 
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EOM scheme has been reported a long time agOjiS a i2L2£ 
but was ignored in more modern uses of it^ which are 
designed to study the Kondo peak in the density of states. 

The failure of the EOM method to describe faithfully 
the Anderson model at its symmetric point, where the 
Fermi level lies exactly between the states of single and 
double occupancies, is a very severe drawback of this 
method. An important question is whether this point 
is singular, or is there a continuous domain in which the 
EOM method fails totally. We return to this problem in 
the next subsection. 



B. Zero-temperature relations 

The zero-temperature limit is of special importance 
since the Green function at the Fermi energy at T = 
satisfies the Fermi-liquid relations 4-29 



ImlG+tO)]- 1 ^, 
RelG+tO)]- 1 = -r CT cot(7rn dCT ) 



(16) 
(17) 



Here and below we use A ± (uj) = lim I) _ = |-o + iv): so 
that G+(w) = [G~(<jj)\* is the usual retarded Green func- 
tion. In Eqs. (|To)) and (|T7I) . T a is the level broadening, 



r (T = -ImS+(0) 



(18) 



The first relation, Eq. (|16|) . implies^ 9 - number conserva- 
tion, and therefore is sometimes referred to as the 'uni- 
tarity' condition. The second one, Eq. (|17fl . is the Friedel 
sum-rule, in which fida is the total number of spin a elec- 
trons introduced by the quantum dotfS 9 - 



nda 



ilm J /(«) 



G+(uj)duj. (19) 



When the self-energy T, a does not depend on the energy, 
fida coincides with the single-spin occupation number on 
the dot, (n da ). 

It is evident that the EOM solution for the Green func- 
tion at the particle-hole symmetric point, Eq. I|15|) . 
olates the unitarity condition Hit) II . On the Fermi level, 
the particle-hole symmetric Green function is 



U 



4(r CT + 2r s 



(20) 



and therefore the imaginary part of [G+(0)] _1 is not de- 
termined solely by the non-interacting self-energy, (as 
implied by the unitarity condition), but has also a con- 
tribution coming from the interaction. This is contrary 
to the result of Lacroixfi 2 - whose EOM differs from our 
Eq. in several places. On the other hand, the Friedel 
sum- rule is satisfied by the Green function l|15l) . which 
yields h d(J = 1/2. This follows from Eq. (Q]|J|: The 
imaginary parts of both G a (ui) and £ CT (u;) are even in u/, 
while (at the symmetric point) the real parts are odd in 
it. However, G a {u>) ~ w _1 at large frequencies, whereas 



dT, a (oj) / dui ~ u~ 2 . As a result, the second term in the 
square brackets of Eq. (|19f) does not contribute. With 
hda = 1/2, the Friedel sum-rule gives Re[G+(0)] _1 = 0, 
which is fulfilled by Eq. (|2T)jl . 

It is rather intricate to study the full EOM solution, 
Eq. at T = 0, even on the Fermi level. However, 
there are cases in which this can be accomplished without 
constructing the full solution. The investigation of these 
cases will also allow us to examine the behavior of G a (0) 
as the particle-hole symmetric point is approached. To 
this end we note that at T — the functions P(oj) and 
Q(lu) acquire logarithmic singularities as lu — > 0, 

p, T ± H~--r CT Gj(o)inM + 0(i), 

Qt(cj) ~ -ir ff [1 + (0)GJ(0)] In M + 0(1) . (21) 

Therefore, we may examine special points at which the 
functions P and Q are divergent, keeping only the diver- 
gent terms in Eq. 10. Then, that equation reduces to 
an algebraic one, which can be easily solved to yield G a 
at those special points. 

Let us first consider the case in which the Zeeman field 
h on the dot vanishes, but 2eo + U ^ 0. Using Eqs. f2*T|l 
in Eq. © yields 

[G+tO)]- 1 + E+(0) = [Gr(O)]- 1 + Sr(0) . (22) 

By writing the Green function in the general form 

[G^z)]- 1 =z-e -ah-Z a (z)- E™*(*) , (23) 

in which E lnt is the self-energy due to the interaction, 
Eq. takes the form 



Sf + (0)-Ef-(0) = 0. 



(24) 



When the network is not spin-polarized, the spin in- 
dices a and a are indistinguishable. Then Eq. (|22|l im- 
plies that ImS lnt (0) vanishes, namely the unitarity con- 
dition is satisfied. (In the more general case of possi- 
bly ferromagnetic leads, it is only the imaginary part of 
[G+W]- 1 - [Gg (0)]- 1 which is determined by the non- 
interacting self-energy alone.) Had we now sent 2eo + U 
to zero, we would have found that the EOM result at 
the particle-hole symmetric point does satisfy the unitar- 
ity condition, in contradiction to our finding, Eq. 12U|) . 
above. We thus conclude that the failure of the EOM 
to obey the Fermi-liquid relation l|16|) at the symmetric 
point is confined to the symmetric point alone, namely, 
the imaginary part of S lnt on the Fermi level has a dis- 
continuity. 

Next we consider the case where 2eo + U = 0, but 
h 7^ 0. Using Eqs. (JjHJ in Eq. © we now find the 
relation 

[G+m- 1 + S+(0) = -{IG+iO)}- 1 + S+(0)} . (25) 

Inserting here expression J23J), we re- write this relation 
in the form 



Tint + / 



(0) + £f + (0) =U. 



(26) 
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Therefore, lm£ lnt (0) = 0, in agreement with the unitar- 
ity condition. Sending now the Zeeman field on the dot 
to zero, yields the result 2Re£ mt (0) = U, which agrees 
with the real part of Eq. l|20|). Thus the EOM result for 
the real part of S lnt on the Fermi energy does not have 
a discontinuity. We hence conclude that the EOM tech- 
niques failure at the symmetric point is confined to the 
imaginary part of the interacting self-energy alone and 
to the symmetric point alone. 

Our considerations in this subsection are confined to 
oj = 0, and therefore do not allow us to investigate the 
Friedel sum-rule easily. We carry out such an analysis 
for the infinite-?/ case below. Alternatively, one may at- 
tempt, as has been done in Ref. EH to impose the Friedel 
sum-rule on the EOM result, assuming that Eq. I|17(l 
holds, with fida = {ndo)- This is a dangerous proce- 
dure, which leads in some cases to un-physical results, as 
is demonstrated in Sec. ^ 

IV. EXACT SOLUTION IN THE U -> oo LIMIT 

In this section we present an exact solution of the self- 
consistently truncated EOM, and obtain the dot Green 
function, in the limit U — > oo. For this solution, we as- 
sume that the bandwidth D a is larger than the other en- 
ergies in the problem (except U). In the next section we 
use this function to calculate several physical quantities, 
and compare the results with those of the Bethe anstaz 
technique and other calculations. For simplicity, we also 
assume that the non-interacting self-energy may be ap- 
proximated by an energy-independent resonance width, 
i.e., 

- T*r ff , (27) 

for all the energies in the band, — D a < u> < D a (The 
extension to the case where there is also an energy- 
independent real part to the self-energy is straightfor- 
ward). This assumption is certainly reasonable for a 
range of energies near the center of the band. However, 
using it for the whole band introduces corrections of or- 
der \id/D a \, thus restricting the solution to \uj/D a \ <C 1. 
Then, Eq. © for the Green function, together with the 
definitions (fTOfl . (fTTfl . and (|T2)) . takes the form 

_1 G±(w + ah) = l- {n ds ) - Pt{u + ah),(28) 



where 

Qf(w) = [Lu-e ±iT a -I±{Lu + ah) 

T *(r ff + T a )P±{u + ah)Y 1 . (29) 

Note the shift of energies by ah, compared to Eq. 
The function I a introduced here contains the Kondo 
singularity^ 

itn = ^ f D '*s ! W) 7 . (so) 

To order ©(w/D^), one has 

where is the digamma function. Equation (|28|l for 
the Green function also contains the function P a , which, 
using the assumption l|27|) . is given by 

p^)=^r^^MM^. (32) 

Physically, all the integrals which contain S CT must be 
calculated between —D a and D a , and the resulting Green 
function is calculated only for energies inside the band, 
\u>\ < D a . However, the integral P^(ui) converges even 
when one takes the limit D a — > oo, because G a (uj) ~ 1/lu 
at large |w| [see e.g. Eq. (J2HJ]- If D a is sufficiently large, 
so that this asymptotic behavior becomes accurate and 
since f{io) ~ 1 for u> < —D a , it is convenient to extend 
the range of this integral (and all the related integrals 
below, unless otherwise specified) to the range — oo < 
oj < oo. This introduces errors of order Y a jD a or u>/ D a 
in the results, which we neglect. 

We have thus found that the equation for involves 
an integral containing the function Gg , and thus the two 
functions Gj and Gg are coupled. In addition, the oc- 
cupations (rider) have to be determined self-consistently 
from the Green functions themselves. Our solution for 
the Green function follows the method introduced in 
Refs. ^3 an d CHI This method allows one to turn the 
integral equations into algebraic ones, at the cost of 
additional quantities which have to be determined self- 
consistently from the Green function. First, one intro- 
duces the functions 



= z-eo + iT*- h{z + ah) - i{T & + T„) ^ / duf /(< f j? g , 

7T J z + ah — lu' 

$ CT (z) = z - e - iT s - I a (z + ah) +i(r s +T a )^ [ cLu' 1 ^'^^ . (33) 

7T J z + an — u> 

Note that ^{u>) is identical to fef(w)] 1 and $>~(uj) = \Q^(u>)] , while &~(u>) and 5>+(w) are different from 
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\Q a (w)] 1 and [Cf^(aj)] 1 , respectively. The knowledge of is sufficient to determine G s , since 

$+ (w) - $- (w) = 2^^ ff /(w + aft) [f + i (r ff + Gr ( W + aft)] . (34) 
Similarly, the functions <&J determine G+, through the relation 

- = 2zr CT /( W + aft) [l - i{T s + T a )G+{uo + ah)] . (35) 

Returning now to Eq. I|28|) . we eliminate the Green functions by using Eqs. (|34|l and l|35|l . and the functions P by 
using the definitions (|33|l . In this way we find 

=X», ^M #feM = X-(,), (36) 
2zl o-/(w + aft) 2iL s j(uj + ah) 

where 

X CT (z) = -i(r CT + IV) (1 - (n^)) + z - e + TV - I s (z + aft) , 

X a (z) = i{T a + T s ){l-(n da )) + z-e -iT s -I a (z + ah). (37) 
So far, we have not achieved much simplification over the original problem at hand. However, noting that 

X+(w)-X-(w) = 2iT s f(Lo + ah) , X+(u)-X~(u) = 2iT a f (uj + ah) , (38) 
Eqs. 1)36(1 yield the remarkable result 

$+ (w) *+ (w) - X+ (w)X+(w) = (w) * - (w) - X; H X; ( w ) . (39) 
Therefore, the combination 

i?(z) = $ CT (z)^(z)-X CT (z)X ff (z) (40) 

is non-singular across the real axis. In fact, the only singular point of this combination is at z = oo. This means that 
R(z) can be written as a polynomial with non-negative powers of z. Moreover, since $, $, X, and X grow only linearly 
as z — > oo, that polynomial includes only two terms, ro + r\z. The details of this calculation are given in Appendix 
IbI The result 139|l allows one to express the (unknown) functions $ and $ in terms of the (known) functions X, X 
and R, 



$+ .R + X+X" , $+ R + X+X+ * , . 

= = Ha(w) , ^ = " = H CT (w) . 41 

*<r i? + X^X^ <S>„ R + X+X^ 

This reduces our problem into two independent linear Reimann-Hilbert problems, for which a rigorous solution is 
available! 31 ! 32 



$ a (z) = (z-a)e M °(*\ M a (z)= f (-^ 

J \ 2m 



duj \ lnH (T (a;) 



= (z-oJe^W, M,W=/(-^)^. (42) 

/ V 2m J z — u> 



This is a valid solution as long as lnH(u;) and lnH(w) can 
be chosen to be continuous in to and to vanish at both 
ends of the integration intervali 3 ^ All the cases studied 
in this paper obey this requirement, although we could 
not prove the absence of solutions other than l|42|l for a 
general case with no spin symmetry. The form of the 
polynomial prefactors, (z — a) and (z — a), in Eq. ((4*2^) is 
dictated by the fact that the leading term in $ a {z) and 



<&cr(z) must be z [see Eqs. (|33|) ]. The determination of 
the coefficients a and a, as well as other self-consistent 
quantities, is detailed in Appendix IU1 
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FIG. 1: The zero-temperature occupation number no as func- 
tion of the renormalized energy Ed calculated self-consistently 
by the EOM method (solid line), and from the Friedel sum 
rule, Eq. 11171 (dashed line). Open circles show the exact 
Bethe ansatz results?— 



V. PHYSICAL PROPERTIES IN THE U -> oo 
LIMIT 




FIG. 2: The local spin susceptibility \ as function of the 
renormalized energy level Ed, for T = 0. The EOM re- 
sult (solid line) is close to the Bethe ansatz one— (circles), 
except for large negative values of Ed/T. Inset: The scal- 
ing of the susceptibility with T/T^ at three fixed energy 
level positions (marked by arrows in the main graph). The 
solid lines indicate the high-temperature asymptotic behav- 
ior, X (T) = ( 9f i B ) 2 /(6T). 



Once the functions $ and $ are found, then the Green 
function is determined from Eq. (|34[) or Eq. (|35[1 . This 
knowledge enables us to compute various physical quan- 
tities, and compare them with the results of other cal- 
culations. The first quantity we consider is the total 
occupation on the dot, no = (n^ at zero tem- 

perature. This calculation is carried out for the spin- 
symmetric case, h = and spin-independent self-energy. 
We also denote D = D a , T = T a . The result is plotted 
as function of Ed/T, where 

E d ee e + (r/7r) HD/T) , (43) 

and is portrayed in Fig. ^ (full line). It agrees within 3% 
with the exact universal curve no(Ed/T), as found from 
the Bethe ansatz^ (open circles). Thus, the EOM so- 
lution conforms with Haldane's scaling^ On the other 
hand, the EOM solution fails to satisfy the Friedel sum- 
rule, as has been already discussed above. The total oc- 
cupation calculated from Eq. (|17fl (dashed line) deviates 
systematically from the self-consistent values, in partic- 
ular in the Kondo regime (no —►!.). Note that Eq. I)27|) 
implies that hd a — (nda) [see Eq. |[TSJ]. 

As mentioned above, the Fermi-liquid relations are 
connected with unitarity. In particular, at zero tempera- 
ture the linear conductance of a symmetrically-coupled 
quantum dot is given by2£ — 2(e 2 /h)T a lmG a (u) = 
0). Thus, the unitary limit 2e 2 /h is reached only if 
Re[G+(0)] — > and the Fermi liquid relation holds. 
As implied by the dashed line in the figure, the first of 
these criteria is not obeyed by the self-consistent solution. 

Next we consider the local spin susceptibility on the 
dot. When the leads are non-magnetic, this quantity 
is given by \ = \(g^ B ) 2 d(n d i - n d ]) / dh, with h = 



—gfisB. Here B is the external magnetic field and g^B 
is the gyromagnetic ratio of electrons in the quantum 
dot. (The spin susceptibility of the leads adds to the lo- 
cal susceptibility the usual Pauli term, and 0{T a /D a ) 
corrections^) We have calculated x(Ed/T,T) by differ- 
entiating the self-consistent equations for {rider) with re- 
spect to h, and evaluating the integrals numerically. This 
procedure is similar to the one which has been used for 
the Wolff model in Ref. II3b, but is free from numerical 
accuracy problems reported there. 

The zero-temperature susceptibility derived from the 
EOM is found to be in a good quantitative agreement 
with the Bethe ansatz results in the mixed valence 
(|i?d/r| < 1) and empty orbital (Ed ^> T) regimes, as is 
shown in Fig. El I n the local moment regime, Ed <C — T, 
a screening cloud is expected to form due to the Kondo 
effect T <T K ~ re 7rBd/(2r) , leading to a crossover 

from a high-temperature Curie law, x ~ (qV-b) 2 /T , to 
a finite ground state value, x ~ (<? u b) 2 /Tk- The latter 
is underestimated by our self-consistent solution, as is 
manifested by the deviation from the exact Bethe ansatz 
results depicted in Fig. [21 

Indeed, had we defined the Kondo temperature 
through the inverse of the zero-temperature susceptibil- 
ity, we would have found that the EOM method overesti- 
mates that temperature. However, within EOM, the rel- 
evant energy scale is determined from the leading (real) 
terms in the denominator of the Green function, i.e., by 
the temperature at which the real part of Q, Eq. Q2'jp. 
vanishes. Using for I(z), Eq. I|31|l . the approximate 
forms I(z) m —(T/tt) log [(z + inT) / D] where k is a num- 



ber of order unity, we find that the leading terms are 

_ T d r v^ 2 + « 2 r 2 

u - E d + - In - + - In , (44) 



yielding for the temperature scale 
T* K ~ Te nEd/T , 
such that the leading terms are 



+ ^in^(o;/r^ + («r/r^ 



(45) 



(46) 



The logarithm in Eq. (|46[) dominates the properties of 
the solution close to the Fermi energy at temperatures 
T < Tjr. The same energy scale T£ has been determined 
from the analysis of the truncated EOM in Refs. and 
E Note that T* K is smaller than the true Kondo tem- 
perature Tk- 

The local spin susceptibility at finite temperatures, cal- 
culated from our EOM solution, is shown in the inset of 
Fig. [2 Indeed, x(T)/x(0) scales with T/T£, but instead 
of crossing over to the Curie law, a region of intermedi- 
ate behavior in which x ~ T x with — 1 < x < is ob- 
served. The high-temperature asymptotic x(T) ~ l/T is 
approached only for T >T. Such a behavior in the inter- 
mediate temperature range T£ < T < T is not supported 
by NRG or Bethe ansatz calculations, 4 which scale with 
Tur^ and has thus to be attributed to the deficiency of 
the EOM method. 

In contrast, neither the Lacroix approximation, as im- 
plemented in Ref . flU nor the MWU approximate Green 
function leads to comparable results when used to calcu- 
late the local spin susceptibility. The Lacroix approxi- 
mation becomes intrinsically inconsistent at finite mag- 
netic fields, since it results in a logarithmic divergence 
of G a (z) as z — ► ah. Even when ignoring this inconsis- 
tency, the zero-temperature spin susceptibility calculated 
in that approximation attains negative and divergent val- 
ues regardless of the quantum dot parameters. The MWL 
approximation leads to finite, but quite unphysical values 
of x for T < r, as we demonstrate in Fig. [21 That approx- 
imation gives reasonable results at T = T: the suscepti- 
bility follows roughly the Curie law « 1/T(= l/T) or the 
zero-temperature value, whichever is smaller. At lower 
temperatures one would have expected a gradual increase 
in the susceptibility in the Kondo region. Instead, a win- 
dow of a negative susceptibility opens, which is widened 
as the temperature is decreased. At strictly zero tem- 
perature x is negative for all values of Ed- This example 
shows the necessity of using the full self-consistency of 
the EOM solution in order to obtain the qualitatively 
correct behavior. 

Finally, we examine the local density of states on the 
dot, given by 



^,T)E-Im^G+( W )/x, 



(47) 




5 e / r 1 ° 

d 

FIG. 3: Data from Fig. |3] compared to the spin susceptibility 
given by the MWL approximation in the temperature range 
from T = ir down to T = 1(T 4 F The inset shows the 
MWL susceptibilities only, on a linear scale, allowing for the 
negative values. 
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as calculated from the EOM. The inset of Fig. 21 shows 
the Kondo peak at temperatures T < Tf ( . The local 



FIG. 4: (Color online) The scaling of the local density of 
states at the Fermi level as function of the reduced tempera- 
ture T/Tk, is demonstrated by plotting p(co — 0,T) for three 
different values of the renormalized single-electron energy on 
the dot, E d [Eq. g3J]. The dashed line at p = 2/(tiT) 
corresponds to the limit dictated by the Friedel sum rule at 
Ed/T <g; — 1. The inset shows the melting of the Kondo peak 
as the temperature is increased at fixed E^/T = —5. Temper- 
ature values represented in the inset are marked by arrows in 
the main graph. 



density of states at the Fermi energy follows the univer- 
sal temperature dependence, as can be seen in Fig. 01 
with the same scaling factor T£ as the spin susceptibil- 
ity (compare to the inset of Fig.[5J). The appearance of a 
single scale determining the low energy properties of the 
system is another hallmark^ of the Kondo effect which is 
captured by the fully self-consistent EOM technique. 
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VI. PERTURBATION EXPANSION IN THE 
DOT-NETWORK COUPLINGS 



where we have denoted 



The EOM technique is unfortunately not a systematic 
expansion. It is therefore very interesting to compare 
its results with those given by a direct expansion. Here 
we expand the EOM Green function up to second-order 
in the dot-network couplings J na . For compactness, we 
confine ourselves to the case U — > oo, and assume for 
simplicity that the non-interacting self-energy has just 
an energy-independent imaginary part, i.e., Y, a (oj) = 
—iT a . By comparing with the direct perturbation the- 
ory expansionist we find that, up to second-order in 
J na , our EOM result for the dot Green function is exact. 
On the other hand, the Green function derived in Ref. 
(which does not include all correlations resulting from the 
truncated EOM) violates the second-order perturbation 
theory result, and in fact, predicts a Kondo anomaly in 
this order (which should not be there). 

The U — ► oo Green function of the EOM is given 
by Eqs. (3$, {SJJ, J23), and Since the non- 

interacting self-energy E CT , and consequently the width 
T a are second-order in the coupling J na , the expansion 
of the Green function reads 



x (°) 



Z - e 



(48) 



(n s ) = 1 - Sn s , e a = e + ah 



(50) 



and where the superscript (k) denotes the contribution 
of order k in the couplings. The function P^ is found 
by using the zeroth-order of the Green function in Eq. 



P^(z) = G^{z)[l a {z)-I-{e a ) 

- /(e CT )S CT (z) + /( £CT )S;(e a )] 



(51) 



Using the identities 



I-(e a ) - /(e CT )E-(e CT ) = Rc[l a (e a ) - f(e a )E„(e a j\ (52) 



and 



and 



z - e a 



Gf (z 1 ) = 5n^G^(z)/84 



(0) 



(49) Eq. g5J takes the form 



(53) 



G&(z 



X I 2 



5n^ ] \h{ Zl ) - /( eg )S g (z 1 ) - Re{/ g (e g ) - /(e g )E g ( £g )}] + 5n^ \I s ( Zl ) + S g ( 

(z - e a f 



(o) 



(54) 



r 



In order to complete the second-order calculation, we 
need to find the occupation numbers Sn^ and drffl [see 
Eq. l|5UJ)]. From Eq. gSJl we find 



(0) 



/(e g )[l-/( eg )] 

1 - /(O/fer) 



(55) 



The second-order correction to the occupation number, 



5ri 



(2) 



\(2) 



is obtained by integrating over the 



Fermi function multiplied by the second-order correction 

(2) 

to the density-of-states, p a (w). The latter reads 



G+ (2) H-G" lz » 
-2ni 



(2), 



<5( W -e g )(fe g 2; + (n„>W-(n ff ) 



v(0) 



\(0) 



dRe/ g (e g ; 
<9e 



^(w-e g )fe g 0) Re/ g (e g 



T ff + (r g - IV) (n g )W + r g /(q; - e CT + e g ) - (n g ) 

tt(uj - e a ) 2 



(0)1 



(56) 
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Apart from the terms representing the second-order modifications of the singularity at e CT [the first and second 
members of Eq. (|56[) ]. our result reproduces the one of Ref. 122. for the case where the width r a is spin- independent, 

i.e., r CT = T s . Note that this density of states remains finite at all temperatures. The Kondo divergence of pt^O) in 
the limit T — > appears only at k > 4 [see Ref. l22j . Using the second-order correction to the density of states, we 
find 



(n.) (2) - J f{u)p£\u)du = ((n ff ) (0) ~ l) Re^^ + ^—Rel^e,) + ^—Rel^) . (57) 




FIG. 5: The equilibrium occupation number no — X^a( na ") 
for the same parameters as in Fig. Q and D = 100r (since 
second-order perturbation theory does not scale with Ed, 
the bandwidth has to be specified), calculated from the self- 
consistently truncated EOM (solid lines) , and by perturbation 
theory to first-order in V, Eqs. 1551 and 1571 (dashed line). 
The perturbational result diverges at eo = 0. 

This result is identical to the U — > oo limit of Eq. (5) 
in Ref. which was obtained by a direct perturba- 
tion expansion. Figure |3] depicts the total occupation, 
no = J2a( n <?)' as function of eo, a s found from the EOM 
technique, and as computed from Eqs. (|55|l and 157|) to 
first order in the width. The two curves differ by a few 
percents. The comparison with the exact result is carried 
out in Sec.lTVI 

VII. CONCLUSIONS 

We have presented a solution for the Green function 
of an interacting quantum dot embedded in a general 
non-interacting network. Our solution is derived within 
the EOM technique, taking into account exactly all re- 
sulting correlations once those equations are truncated. 
We have tested our solution by analyzing several limit- 
ing cases, and by comparing several physical properties 



I 

derived from that solution with other results available by 
the Bethe ansatz method and by NRG computations. 

We have found that the EOM Green function is 
temperature-independent at the particle-hole symmetric 
point (where h = and 2eo + U = 0). We have found 
that this deficiency of the EOM is related to a disconti- 
nuity in the imaginary part of the interacting self-energy 
at that particular point. However, this imaginary part 
obeys the Fermi liquid unitarity requirement away from 
this special point, at zero temperature. In contrast, even 
away from the particle-hole symmetric point, the EOM 
result fails to satisfy the Friedel sum-rule deep inside the 
Kondo regime, as we have shown explicitly in the infinite 
interaction limit. 

Albeit these problems, the EOM solution reproduces 
faithfully the low-temperature scaling of the spin sus- 
ceptibility and the density of states at the Fermi level, 
though with an energy scale Tj*- which differs from 
the true Kondo temperature [ cf. Eq. I|45fl]. Zero- 
temperature results are in excellent agreement with the 
exact Bethe ansatz solution, except for the Kondo corre- 
lated regime. As the temperature is raised, the EOM 
results become more and more quantitatively correct, 
and approach high-temperature asymptotics known rig- 
orously from perturbation theory and NRG studies. We 
have expanded the EOM Green function to second order 
in the dot-network coupling and found an exact agree- 
ment with direct calculations by perturbation theory. 
Most importantly, we have found that it is crucial to in- 
clude in the EOM solution all the correlations emerging 
from the truncated scheme. Ignoring part of these corre- 
lations, as is ubiquitously done in such studies, results in 
erroneous behaviors of various physical quantities. 

Hence we conclude that the method examined in this 
paper can provide a reasonable description of a quantum 
dot system over a wide parameter range, provided that 
the self-consistency conditions inherent to this technique 
are fully taken into account. 

We thank A. Schiller for helpful comments. This 
project was carried out in a center of excellence supported 
by the Israel Science Foundation. 
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APPENDIX A: EOM FOR THE FINITE U CASE 

Here we extend the derivations of Rcfs. HE1 andE carried out for an infinite repulsive interaction, to the case in 
which U is finite. In addition, we allow for a Zeeman field, assuming that the quantization axes of the spins are the 
same on the dot and on the leads (but the g-factors may be different). The Hamiltonian of the our system is given in 
Eq. CJ. 

The Fourier transform of the EOM for the Green function defined in Eq. JSJ) can be written in two alternative 
formal (both will be used in the following) 

z<A;B> 2 = ([A,B] + ) + <£[A,H]_;B^> z = ([A,B] + ) - «A; [B,H]_^> Z . (Al) 

It follows that the EOM for the dot Green function is 

[z-e - ah] G a = 1 + U<^n ds d G \ 4> - Y J n^a n(T ; 4> , (A2) 

71 

where a is the spin direction opposite to a. Here and in the following we frequently omit for brevity explicit indications 
of the argument z. We first consider the last term on the right-hand-side of Eq. (|A2|) . The EOM for the Green function 
appearing there is 

(z - e na ) <a Ilcr ; 4> = - Y ^nm^ma! 4> - Jn*Ga ■ (A3) 

m 

Introducing the inverse matrix 

M nmcr {z)=[{z-Hl e% )-\ m , (A4) 
where 7i^ ct is the part of Ti nc t pertaining to the spin direction <r, we find 

«a„ ff ;4» = -S^ n Jz)CG (F . (A5) 

m 

Note that M. is the Green function matrix of the network in the absence of the coupling to the dot. Using Eq. i|A5|l , 
we find that the last term on the right-hand-side of Eq. I|A2() can be put in the form 



- Y J n<r<-a n(7 ; 4> = Z<? G * » ( A 6) 

n 

where 

£ CT (z) = Y J n «M nma {z)j; n!J (A7) 

nm 

is the self-energy of the dot Green function coming from the coupling to the (non-interacting) network. Namely, it is 
the dot self-energy for the U = case^ As such, it can always be calculated, at least in principle (see for example 
Refs.GJandizI). 

We now turn to the interacting part of the EOM for the dot Green function [the second term on the right-hand-side 
of Eq. ilA2(l ] . The EOM for the 4-operator Green function appearing there reads 

[z-e -ah- U] <n da d a ; 4> = ( n ds) - Y [ J na<a n<T n da ; 4> 

71 

+ J nff <4a„a4; 4> - <C<aL?44; 4> - ( A «) 

and gives rise to three new 4-operators Green functions (on the right-hand-side here). Their EOM are 

[z - e na \ <a„ CT n dff ; 4> = - Y J ™n<.a mrj n ds \ 4> - Ka<- n ds d a\ 4> 

m 

- Y [ J ™ff <a mr4a m( r;4> ~ <Ca< a Lx4 a mx 5 4> . ( A9 ) 
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[zi - e nS ] <^dla nS d a ; 4> = (4*0 - J3 J nm <4<w4; 4> - Jns<- n d S d a ] 4> 

m 

+ X] [ J mff< a La- a na4; 4> ~ 4 ff «4a nf a mff ;4» . ( Al °) 
m 

[-z 2 + e„ ff ] <4 5 d 5 d CT ; 4> = (4ffO + X! J nm< a ma44;4> + Jna<n d ^d a ; 4> 

m 

- ^ [j TOCT <4 s C^a TO(7 ;4>- ^ma<4a4 a TOa;4> > ( A11 ) 
m 

where we have introduced the definitions 

zi = z - 2cr/i, z 2 = -z + 2e + C7 . (A12) 

The EOM l|X9)l - l|ATT|l include 4-operator Gr ccn functions in which only two of the operators are dot operators. Those 
are decoupled as detailed in Eq. JSJ. One then finds 

Jn<T<an<r"5;4> =E <T <n s d .;4>+ ^ J rwMnm<r<a m<7 ;4> 4'}(4Vj) _ Jm>A a m'adg) , (A13) 
n nmm' 

- 53 J nff<4a„a4; 4> =(1 + S CT G CT )P ff (21) + S s (z 1 )<n ff d CT ; 4> - G CT Q 5 (zi) , (A14) 

X ^nff<4a44; 4> =(! + S (T G CT )P ff (z 2 ) - Sg.(z 2 )<n s 4; 4> + G a Q^{z 2 ) , (A15) 

where we have introduced 

P (J (z) = - 23 «/«<TA^ m „ CT (z)(4 a m CT ) = - '^{a[ na d c ,)M mn< j(z) J* a , (A16) 

nm mn 

Qo{ z ) = ^ ] JncrM-nma{z){f^ m i a a ma) m' a = ^ ] Jm' a (o-l na a m , a ) A4 mna (z) J na . (A17) 
nmm' nmm' 

The second equality in each of Eqs. I|A16|) and (|A17|I is justified below. 

Examining Eqs. (| A13|> . IjAlfijl and (|A1 T|> reveals that one needs to find thermal averages of two types, the ones 
belonging to two network operators, (aj ricr a m , cr ), and the ones consisting of a dot and a network operator, {d\a ma ). 
These are found [see Eq. J7J] from the corresponding Green's functions, whose EOM are given by Eq. (| A3|) and 

(z - e na ) <d CT ; 4 CT > = - ^ <d CT ; 4 lCT >J m „ - JnaG a , (A18) 

m 

{z - e ma ) <a mCT ; a\ CT > = S mn - ^ 4m'«o m v; a L» ~ J m^ d a^ a L> ■ ( Al9 ) 

m' 

Their solutions in terms of the inverse matrix A4 nma [see Eq. (|A4|) ] are 

<d (T ; 4 ff > = - J3 -/maMmno-G,, , (A20) 

m 

<a mcr ; 4 CT > = Mmna + ^3 ■^mlaJ^ J n' aMn'naG a . (A21) 



in' 



Employing these solutions to obtain the thermal averages appearing in Eq. l|A13Jl . we find Yl m ' Jrn's{d\a m , s ) = 
Ylm' Jm's{ a m 'a d a)- Consequently, the terms in the square brackets of Eq. (|A13|) are cancelled. Next, we use the first 
equality in each of the definitions (|A16(1 and i|A17(l together with the auxiliary Green functions (|A5|) and (|A21(1 . to 
find P and Q in terms the dot Green function G, 

P °^ = i <f ^f( z ') G °( z ')Y, J ™ J *nA(z-K c tr\z-K c tr 1 } nm , (A22) 

JC nm 

Q"(*) =i fJ^M [l + S^zOG.^OlE^^^-^ncJ^^'-Kct)- 1 ]^ , (A23) 
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where we have used Eq. (|A4|) . Since 



(z-nur'iz-ni 



-1 _ i z ' ^nct) 1 ( z ^net) 



(A24) 



we can use Eq. (|A7J| for the non-interacting self-energy, to write the functions P and Q in terms of that self-energy, 



P a {z) = lim — / doj 

ri'—>0 Zir J z — oj 



(g o {oj + irf) [E^w + irf) - S CT (z)] G a (co «/) [E a {uj - irf) - £„(«)]) 



}{w)G <J {w) dw, 

27T ./C z — w 



Q a (z) = lim J- / dw- /( "'' 



z — OJ 



x (jl + £ CT (w + irf)G a {u) + if]')] [T, a (u + ir]') - T> a {z)\ - [1 + S CT (w - irf)G a {oj - it]')] [£ CT (u> - «/) - £<t(.z)]) 

(A25) 



= ^ / f(w)[l + X a (w)G a (w) f* {w) ^ {z) dw. 
2vr Jc z-w 



r 



Here and elsewhere the imaginary part of z is always 
greater than rf, so that the contour C never encircles 
the pole at w = z. Note that using the same procedure 
employing the second equalities in Eqs. i|A16(l and (|A17|) 
gives again Eqs. I|A22|I and l|A23[) . thus proving that the 
two definitions of P and Q in Eqs. I|A16|I and l)A17j) are 
equivalent. Equation 1|A8(I can be now easily solved. In- 
serting the solution into Eq. (|A2|) leads to the expression 
for the dot Green function, Eq. (JHJ). 

In treating the functions P and Q in Sec. IIII Al we have 
employed several properties of the complex integrals ap- 
pearing in Eqs. I|A25|I . Consider for example the integral 



(A26) 



Since G{w) and £(u>) have no singularities except for a 
cut along the real axis, it is expedient to complete each 
half of the contour C by a large-radius semi-circle in the 
upper and lower half-planes. Then, by the residue the- 
orem, Eq. (|A26() vanishes provided that G(w) falls as 
w^ 1 or faster at w — > oo. This means that the Fermi 
function f(w) in the definition of P can be replaced by 
f(w) + const. Another important case is when G(w) is 



J 



replaced by 1 in Eq. (|A26|I . In this case the contribu- 
tion of the semi-circles does not vanish, and the integral 
iHH gives £(z). 

An alternative to the fully self-consistent treatment in- 
vestigated in this paper has been proposed in Ref. 5. 
There, the averages of the form (dJ T a miT ) were ignored, 



and those of the type 



"ma m u i 



were approximated by 



7n(T 711(7 1 



Namely, 7i nc t-dot was put to zero in the 
calculation of the averages. Upon such an approxima- 
tion, P<j{z) ~ 0, and 



Im 



- irf) 



du 



(A27) 



reproducing Eq. (8) of Ref. 0. Another approximation of 
Eqs. (|A22|> and l|A23fl . originally due to Lacroixjii has 
been recently analyzed in the context of quantum dots 
in Ref. 0. In this approximation one assumes the dot 
Green function to vary smoothly enough over the inte- 
gration regimes in Eqs. I|A22|) and (|A23(I . so that it can 
be taken out of the integrals. This ansatz reduces Eq. (JSJ 
of the text to a quadratic form. These two approximate 
solutions are discussed in Sees. IIIII and IIVI 



APPENDIX B: DERIVATION OF THE POLYNOMIAL FUNCTION R 



As explained in Sec. IIVI the function 

R{z) = $ CT (z)$ CT (z)-X CT (z)X CT (z) 

= i(T a + Ta) z — eo — iTa - I a (z + ah) 

- i(T a + Ta) z - e + «r CT - I„{z + ah) 



Y- a f f(uj)G^(uj) 

1 - n ds ) / duj — — — 

7r J z + ah — oj 

r CT r f(oj)Gt(uj) 

1 - n da ) / dw — ■ — 

tt z + ah — oj 



(F CT + Tg-J (1 - n ds ) (1 - n da ) / dw — x — f du 



z + ah — oj it 



z + an — oj 



(Bl) 
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is non-singular across the real axis, and its only singular point is at z = oo. This observation enables one to solve for 
the Green function in terms of the functions $ and $. Here we examine R(z) in some detail, and also derive the first 
two terms of its polynomial expansion. 

In the limit z — > oo, the function I&(z), Eq. is given by 



r d 

TT Z 



(B2) 



The z — > oo limit of the integrals appearing in Eq. (|B1|I has to be taken with care. We write the Green 
functions appearing in the integrands there in the form G^(u>) = ReG+(w) ± ilm G~£(u)), and use {rider) = 
— (I/71") J dtaf(u))lmG a ~(u)) to obtain 



r 

TT 



— / duj 



where 



4r(*) 



7T 



do; 



/(«)HeG+(w) 



(B3) 



(B4) 



Inserting Eqs. (|B2|I and l|B3|l into Eq. (|B 1 1) . the terms which survive the z — > 00 limit are 

i?(z) — (r CT + T s ) 2 ((n da ) + (n da ) - (rider) (rids)) + i(T a + Ts)(z - e ) (n da - rids) 
+ i(r ff + r ff ) (z(A a (z) - A s (z)) - (1 - ^ In ^ + (1 - n da ) ^ In ^ 

V IT Z IT Z 



(B5) 



According to the discussion in Sec. IIVI the terms logarithmic in z have to disappear. The integral giving A a (z), Eq. 
(|B4ll . is well behaved on the positive u> axis, since then for large u> the Fermi function makes it convergent. For very 
large negative u> values, ReG^(c<j) — > (1 — rids) /w, and as a result, the contribution from that part of the integration 
to A (T (z) is (1 — n c iff)(T a -/'!T)(X/ z) ln(C CT /z), where C, a < _D CT . Hence, the terms logarithmic in z are cancelled. In our 
calculations, we have used 



— (1 - ri ds ) ln-^- 

TT D a 



A a {z) 

and have determined the coefficient b a self-consistently (see next Appendix). In this way we find 

R(z) = (T ff + Y s ) 2 (nda) + (rids) - {nda) (rids) + i(T a + T a ) (z - e ) (n da - n ds ) +b a - b s 



(B6) 



(B7) 



APPENDIX C: DETAILS OF THE EXACT SOLUTION 



This Appendix is devoted to the analysis of the exact solution for the self-consistent EOM, and in particular to the 
determination of the unknown coefficients b a [see Eqs. <|B6(1 and (|B7ll ]. a and a [see Eqs. (|42f) ]. This is accomplished 
by expanding the solution at large frequencies, and equating the coefficients with those of the desired functions $ and 
$. We give the details for $; those of the 'tilde' solution are obtained analogously. 

As has been the case for the integral (|B3|) . the large negative part of the integral defining M a (z) has to be taken 
with care. To this end we write 



**>-y(-5i) 



duj \ In Ha(u>) 



duj \ InHaiuj) - Q(-D s - u)Fa{u){-2m) 
liri) z — ui 



Fa{u>) dui 



(CI) 



where the function F a (u}) is defined in such a way that the use of the geometric series l/{z — u>) = z^ 1 + z~ 2 oj + 0(z~ 3 ) 
in the first integral of Eq. IjCljl results in convergent integrals. It is sufficient to include in F a {ui) the most slowly- 
decaying terms of hiH(u>)/(— 2iri), which are obtained by expanding H(u>) for large negative u>, 



Fa(u) 



Ta_ 

ITU) 



— In — y + e - i(T„ + T s )(nda) + «r a 

TT \U!\ 



(C2) 
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[In this expansion one has to include terms of the order 0(u> 2 ) because of the linear term in the prefactor in Eq. 
(I42[) .] Using Eq. (|C2|) . the second integral in Eq. (|C1|) is obtained explicitly, and then expanded up to order z~ 2 , 



-D s 



F a (u)) duj _ r g jz 1 / z 

z — ui irz D s 2 \irz D s 



7T 2 D a Z 6 Z 2 



a / [eo - i(T a + T a )(n dcr ) + iT s ] 

ir Z z 7T 



1 1 



As a result, the asymptotic expansion of the function M a (z) becomes 

M a (z) + (Tff/Tr) \n{z/D s )]/z+ 

{P a + (Tff/Tr) [e - *(r ff + r ff )( ?W ) + «r ff ] ln(z/£> ff ) - (r CT /7r) 2 (l/2) bx 2 (z/D 9 )}/^ , 

where the coefficients a<j and j3 a are given by 

p2 p ;-+oo 

a - = ^ ^-[^-i{T a + T s ){n da )+iT g }+ / [ilni7 ff H/(2^)-e(-^-^)F ff H]da;, 

r? r-D- /"+ 00 

/3 ff = --^ + ^^+ / w[ilniT (T H/(27r)-e(-D s -cj)F (T (w)]da;. 

6 7T J_ oc 

An analogous calculation gives F a (uj) = Fg.(w)*, leading to 

M a {z) + (T^/tt) ln^/A,)]/^ 

{& + (r CT /vr) [e + i(T a + T s )(n ds ) - iT a ] \n{z/D a ) - (TJtt) 2 (1/2) ln 2 (z/D a )}/z 2 , 

r 2 r„ 



7t 2 D a IT D„ 



[e + i(r ff + r ff )(n dff ) - iT a ] + 



ilnfl" CT (w)/(27r) - 0( D a - ui)F a (w) 



dio , 



6 



i hiH^u)/ (27r) — 6 (—Do- - u)F a {uj) 



dtu . 



(C3) 



(C4) 

(C5) 
(C6) 



(C7) 
(C8) 

(C9) 



Finally, we use Eqs. I|C4|) and l|C7() in Eq. (|42ll . and then compare term by term with the expansion of Eq. J22J>- 
This procedure determines the coefficients a and a, 



a = a a + e - iT a , a = a a + e + lT CT , 
and gives the self-consistency equations 

b s + iT s {n ds ) = i(T a + Tg)^ 1 [(3 a - aa a + a 2 /2 + T a h/(2%)] , 
b a - ir<j(n dc ,) = —i (r CT + r CT ) _1 [/? CT -aa a + a 2 /2 - r a /i/(27r)] 

I 



(CIO) 

(Cll) 
(C12) 



In the case of a full spin-symmetry (including h = 0), 
one has 5 = a *, ft = ft*, and a = a* , and then Eq. (|C11|) 
and Eq. (|C12(I become complex conjugate. 

For a given set of parameters, Eqs. IjClljl and 1C12|) 



are solved numerically for 6 CT and {n da ) by an iterative 
Newton-Raphson algorithm. The initial values are cho- 
sen from the solution of the non-interacting (U = 0) 
problem. 
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